Localized Soft Vibrational Modes and Coherent Structural Phase Transformations in Rutile TiO2 Nanoparticles under Negative Pressure.

We study the effect of size on the vibrational modes and frequencies of nanoparticles, by applying a newly developed, robust, and efficient first-principles-based method that we present in outline. We focus on rutile TiO2, a technologically important material whose bulk exhibits a softening of a transverse acoustic mode close to q=(12,12,14), which becomes unstable with the application of negative pressure. We demonstrate that, under these conditions, nanoparticles above a critical size exhibit unstable localized modes and we calculate their characteristic localization length and decomposition with respect to bulk phonons. We propose that such localized soft modes could initiate coherent structural phase transformations in small nanoparticles above a critical size.

Experimental signature of a topological quantum dot.

Topological insulator nanoparticles (TINPs) host topologically protected Dirac surface states, just like their bulk counterparts. For TINPs of radius <100 nm, quantum confinement on the surface results in the discretization of the Dirac cone. This system of discrete energy levels is referred to as a topological quantum dot (TQD) with energy level spacing on the order of Terahertz (THz), which is tunable with material-type and particle size. The presence of these discretized energy levels in turn leads to a new electron-mediated phonon-light coupling in the THz range, and the resulting mode can be observed in the absorption cross-section of the TINPs. We present the first experimental evidence of this new quantum phenomenon in Bi2Te3 topological quantum dots, remarkably observed at room temperature.

Correction: Atomistic QM/MM simulations of the strength of covalent interfaces in carbon nanotube-polymer composites.

Correction for 'Atomistic QM/MM simulations of the strength of covalent interfaces in carbon nanotube-polymer composites' by Jacek R. Golebiowski et al., Phys. Chem. Chem. Phys., 2020, 22, 12007-12014, DOI: 10.1039/d0cp01841d.

Atomistic QM/MM simulations of the strength of covalent interfaces in carbon nanotube-polymer composites.

We investigate the failure of carbon-nanotube/polymer composites by using a recently-developed hybrid quantum-mechanical/molecular-mechanical (QM/MM) approach to simulate nanotube pull-out from a cross-linked polyethene matrix. Our study focuses on the strength and failure modes of covalently-bonded nanotube-polymer interfaces based on amine, carbene and carboxyl functional groups and a [2+1] cycloaddition. We find that the choice of the functional group linking the polymer matrix to the nanotube determines the effective strength of the interface, which can be increased by up to 50% (up to the limit dictated by the strength of the polymer backbone itself) by choosing groups with higher interfacial binding energy. We rank the functional groups presented in this work based on the strength of the resulting interface and suggest broad guidelines for the rational design of nanotube functionalisation for nanotube-polymer composites.

The ONETEP linear-scaling density functional theory program.

We present an overview of the onetep program for linear-scaling density functional theory (DFT) calculations with large basis set (plane-wave) accuracy on parallel computers. The DFT energy is computed from the density matrix, which is constructed from spatially localized orbitals we call Non-orthogonal Generalized Wannier Functions (NGWFs), expressed in terms of periodic sinc (psinc) functions. During the calculation, both the density matrix and the NGWFs are optimized with localization constraints. By taking advantage of localization, onetep is able to perform calculations including thousands of atoms with computational effort, which scales linearly with the number or atoms. The code has a large and diverse range of capabilities, explored in this paper, including different boundary conditions, various exchange-correlation functionals (with and without exact exchange), finite electronic temperature methods for metallic systems, methods for strongly correlated systems, molecular dynamics, vibrational calculations, time-dependent DFT, electronic transport, core loss spectroscopy, implicit solvation, quantum mechanical (QM)/molecular mechanical and QM-in-QM embedding, density of states calculations, distributed multipole analysis, and methods for partitioning charges and interactions between fragments. Calculations with onetep provide unique insights into large and complex systems that require an accurate atomic-level description, ranging from biomolecular to chemical, to materials, and to physical problems, as we show with a small selection of illustrative examples. onetep has always aimed to be at the cutting edge of method and software developments, and it serves as a platform for developing new methods of electronic structure simulation. We therefore conclude by describing some of the challenges and directions for its future developments and applications.

Combining Embedded Mean-Field Theory with Linear-Scaling Density-Functional Theory.

We demonstrate the capability of embedded mean-field theory (EMFT) within the linear-scaling density-functional-theory code ONETEP, which enables DFT-in-DFT quantum embedding calculations on systems containing thousands of atoms at a fraction of the cost of a full calculation. We perform simulations on a wide range of systems from molecules to complex nanostructures to demonstrate the performance of our implementation with respect to accuracy and efficiency. This work paves the way for the application of this class of quantum embedding method to large-scale systems that are beyond the reach of existing implementations.

First-Principles Study of Ferroelastic Twins in Halide Perovskites.

We present an ab initio simulation of 90° ferroelastic twins that were recently observed in methylammonium lead iodide. There are two inequivalent types of 90° walls that we calculate to act as either electron or hole sinks, which leads us to propose a mechanism for enhancing charge carrier separation in photovoltaic devices. Despite separating nonpolar domains, we show these walls to have a substantial in-plane polarization of ∼6 µC cm-2, due in part to flexoelectricity. We suggest this in turn could allow for the photoferroic effect and create efficient pathways for photocurrents within the wall.

Multiscale simulations of critical interfacial failure in carbon nanotube-polymer composites.

Computational investigation of interfacial failure in composite materials is challenging because it is inherently multi-scale: the bond-breaking processes that occur at the covalently bonded interface and initiate failure involve quantum mechanical phenomena, yet the mechanisms by which external stresses are transferred through the matrix occur on length and time scales far in excess of anything that can be simulated quantum mechanically. In this work, we demonstrate and validate an adaptive quantum mechanics (QM)/molecular mechanics simulation method that can be used to address these issues and apply it to study critical failure at a covalently bonded carbon nanotube (CNT)-polymer interface. In this hybrid approach, the majority of the system is simulated with a classical forcefield, while areas of particular interest are identified on-the-fly and atomic forces in those regions are updated based on QM calculations. We demonstrate that the hybrid method results are in excellent agreement with fully QM benchmark simulations and offers qualitative insights missing from classical simulations. We use the hybrid approach to show how the chemical structure at the CNT-polymer interface determines its strength, and we propose candidate chemistries to guide further experimental work in this area.

The role of molybdenum in suppressing cold dwell fatigue in titanium alloys.

We test a hypothesis to explain why Ti-6242 is susceptible to cold dwell fatigue (CDF), whereas Ti-6246 is not. The hypothesis is that, in Ti-6246, substitutional Mo-atoms in α-Ti grains trap vacancies, thereby limiting creep relaxation. In Ti-6242, this creep relaxation enhances the loading of grains unfavourably oriented for slip and they subsequently fracture. Using density functional theory to calculate formation and binding energies between Mo-atoms and vacancies, we find no support for the hypothesis. In the light of this result, and experimental observations of the microstructures in these alloys, we agree with the recent suggestion (Qiu et al. 2014 Metall. Mater. Trans. A45, 6075-6087. (doi:10.1007/s11661-014-2541-5)) that Ti-6246 has a much smaller susceptibility to CDF because it has a smaller grain size and a more homogeneous distribution of grain orientations. We propose that the reduction of the susceptibility to CDF of Ti-6242 at temperatures above about 200°C is due to the activation of 〈c+a〉 slip in 'hard' grains, which reduces the loading of grain boundaries.

Chemically Selective Alternatives to Photoferroelectrics for Polarization-Enhanced Photocatalysis: The Untapped Potential of Hybrid Inorganic Nanotubes.

Linear-scaling density functional theory simulation of methylated imogolite nanotubes (NTs) elucidates the interplay between wall-polarization, bands separation, charge-transfer excitation, and tunable electrostatics inside and outside the NT-cavity. The results suggest that integration of polarization-enhanced selective photocatalysis and chemical separation into one overall dipole-free material should be possible. Strategies are proposed to increase the NT polarization for maximally enhanced electron-hole separation.

Unravelling the Roles of Size, Ligands, and Pressure in the Piezochromic Properties of CdS Nanocrystals.

Understanding the effects of pressure-induced deformations on the optoelectronic properties of nanomaterials is important not only from the fundamental point of view but also for potential applications such as stress sensors and electromechanical devices. Here, we describe the novel insights into these piezochromic effects gained from using a linear-scaling density functional theory framework and an electronic enthalpy scheme, which allow us to accurately characterize the electronic structure of CdS nanocrystals with a zincblende-like core of experimentally relevant size. In particular, we focus on unravelling the complex interplay of size and surface (phenyl) ligands with pressure. We show that pressure-induced deformations are not simple isotropic scaling of the original structures and that the change in HOMO-LUMO gap with pressure results from two competing factors: (i) a bulk-like linear increase due to compression, which is offset by (ii) distortions and disorder and, to a lesser extent, orbital hybridization induced by ligands affecting the frontier orbitals. Moreover, we observe that the main peak in the optical absorption spectra is systematically red-shifted or blue-shifted, as pressure is increased up to 5 GPa, depending on the presence or absence of phenyl ligands. These heavily hybridize the frontier orbitals, causing a reduction in overlap and oscillator strength, so that at zero pressure, the lowest energy transition involves deeper hole orbitals than in the case of hydrogen-capped nanocrystals; the application of pressure induces greater delocalization over the whole nanocrystals bringing the frontier hole orbitals into play and resulting in an unexpected red shift for the phenyl-capped nanocrystals, in part caused by distortions. In response to a growing interest in relatively small nanocrystals that can be difficult to accurately characterize with experimental techniques, this work exemplifies the detailed understanding of structure-property relationships under pressure that can be obtained for realistic nanocrystals with state-of-the-art first-principles methods and used for the characterization and design of devices based on these and similar nanomaterials.

Pressure-Induced Amorphization and a New High Density Amorphous Metallic Phase in Matrix-Free Ge Nanoparticles.

Over the last two decades, it has been demonstrated that size effects have significant consequences for the atomic arrangements and phase behavior of matter under extreme pressure. Furthermore, it has been shown that an understanding of how size affects critical pressure-temperature conditions provides vital guidance in the search for materials with novel properties. Here, we report on the remarkable behavior of small (under ~5 nm) matrix-free Ge nanoparticles under hydrostatic compression that is drastically different from both larger nanoparticles and bulk Ge. We discover that the application of pressure drives surface-induced amorphization leading to Ge-Ge bond overcompression and eventually to a polyamorphic semiconductor-to-metal transformation. A combination of spectroscopic techniques together with ab initio simulations were employed to reveal the details of the transformation mechanism into a new high density phase-amorphous metallic Ge.

A first-principles study of the vibrational properties of crystalline tetracene under pressure.

We present a comprehensive study of the hydrostatic pressure dependence of the vibrational properties of tetracene using periodic density-functional theory (DFT) within the local density approximation (LDA). Despite the lack of van der Waals dispersion forces in LDA we find good agreement with experiment and are able to assess the suitability of this approach for simulating conjugated organic molecular crystals. Starting from the reported x-ray structure at ambient pressure and low temperature, optimized structures at ambient pressure and under 280 MPa hydrostatic pressure were obtained and the vibrational properties calculated by the linear response method. We report the complete phonon dispersion relation for tetracene crystal and the Raman and infrared spectra at the centre of the Brillouin zone. The intermolecular modes with low frequencies exhibit high sensitivity to pressure and we report mode-specific Grüneisen parameters as well as an overall Grüneisen parameter [Formula: see text]. Our results suggest that the experimentally reported improvement of the photocurrent under pressure may be ascribed to an increase in intermolecular interactions as also the dielectric tensor.

Simulations of nanocrystals under pressure: combining electronic enthalpy and linear-scaling density-functional theory.

We present an implementation in a linear-scaling density-functional theory code of an electronic enthalpy method, which has been found to be natural and efficient for the ab initio calculation of finite systems under hydrostatic pressure. Based on a definition of the system volume as that enclosed within an electronic density isosurface [M. Cococcioni, F. Mauri, G. Ceder, and N. Marzari, Phys. Rev. Lett. 94, 145501 (2005)], it supports both geometry optimizations and molecular dynamics simulations. We introduce an approach for calibrating the parameters defining the volume in the context of geometry optimizations and discuss their significance. Results in good agreement with simulations using explicit solvents are obtained, validating our approach. Size-dependent pressure-induced structural transformations and variations in the energy gap of hydrogenated silicon nanocrystals are investigated, including one comparable in size to recent experiments. A detailed analysis of the polyamorphic transformations reveals three types of amorphous structures and their persistence on depressurization is assessed.

Bromophenyl functionalization of carbon nanotubes: an ab initio study.

We study the thermodynamics of bromophenyl functionalization of carbon nanotubes with respect to diameter and metallic/insulating character using density-functional theory (DFT). On the one hand, we show that the functionalization of metallic nanotubes is thermodynamically favoured over that of semiconducting ones, in agreement with what binding energy calculations previously suggested. On the other hand, we show that the activation energy for the grafting of a bromophenyl molecule onto a semiconducting zigzag nanotube ranges from 0.72 to 0.75 eV without any clear diameter dependence within numerical accuracy. This implies that this functionalization is not selective with respect to diameter at room temperature, which explains the contradictory experimental selectivities reported in the literature. This contrasts with what is suggested by the clear diameter dependence of the binding energy of a single bromophenyl molecule, which ranges from 1.52 eV for an (8, 0) zigzag nanotube to 0.83 eV for a (20, 0) zigzag nanotube. Also, attaching a single bromophenyl to a nanotube creates states in the gap close to the functionalization site. It therefore becomes energetically favourable for a second bromophenyl to attach close to the first one on semiconducting nanotubes. The para configuration is found to be favoured for resulting bromophenyl pairs and their binding energy is found to decrease with increasing diameter, ranging from 4.35 eV for a (7, 0) nanotube to 2.26 eV for a (29, 0) nanotube. An analytic form for this radius dependence is derived using a tight binding Hamiltonian and first order perturbation theory. The 1/R(2) dependence obtained (where R is the nanotube radius) is verified by our DFT results within numerical accuracy. Finally, bromophenyl pairs are shown to be favoured by only 50 meV with respect to separate moieties on (9, 0) metallic nanotubes, which suggests that pair formation is not significantly favoured on some metallic nanotubes. This result explains the observation of stable isolated moieties at room temperature in nanotube samples containing random nanotube chiralities.

Ab initio calculations of the optical absorption spectra of C60-conjugated polymer hybrids.

A recently developed linear-scaling density-functional theory (LS-DFT) formalism is used to calculate optical absorption spectra of hybrids of C60 and the conjugated polymers poly(para-phenylene) (PPP) and poly(para-phenylene vinylene) (PPV). The use of a LS formalism allows calculations on large systems with realistic proportions of C60, which has been of interest for the use of such materials in photovoltaics. Two different bonding structures are tested for the hybrid PPP and for both systems additional peaks are present in the absorption spectra below the original onset of absorption. By identifying the eigenstates involved in the relevant transitions, a weighted density difference is formed, demonstrating the transfer of charge between the polymer chain and the C60, in agreement with experiment. For the hybrid PPV, no additional peaks are observed in the absorption spectrum.

Electrostatic considerations affecting the calculated HOMO-LUMO gap in protein molecules.

A detailed study of energy differences between the highest occupied and lowest unoccupied molecular orbitals (HOMO-LUMO gaps) in protein systems and water clusters is presented. Recent work questioning the applicability of Kohn-Sham density-functional theory to proteins and large water clusters (Rudberg 2012 J. Phys.: Condens. Matter 24 072202) has demonstrated vanishing HOMO-LUMO gaps for these systems, which is generally attributed to the treatment of exchange in the functional used. The present work shows that the vanishing gap is, in fact, an electrostatic artefact of the method used to prepare the system. Practical solutions for ensuring the gap is maintained when the system size is increased are demonstrated. This work has important implications for the use of large-scale density-functional theory in biomolecular systems, particularly in the simulation of photoemission, optical absorption and electronic transport, all of which depend critically on differences between energies of molecular orbitals.

Toward Ab Initio Optical Spectroscopy of the Fenna-Matthews-Olson Complex.

We present progress toward a first-principles parametrization of the Hamiltonian of the Fenna-Matthews-Olson pigment-protein complex, a molecule that has become key to understanding the role of quantum dynamics in photosynthetic exciton energy transfer. To this end, we have performed fully quantum mechanical calculations on each of the seven bacteriochlorophyll pigments that make up the complex, including a significant proportion of their protein environment (more than 2000 atoms), using linear-scaling density functional theory exploiting a recent development for the computation of excited states. Local pigment transition energies and interpigment coupling between optical transitions have been calculated and are in good agreement with the literature consensus. Comparisons between simulated and experimental optical spectra point toward future work that may help to elucidate important design principles in these nanoscale devices.

Electrostatic interactions in finite systems treated with periodic boundary conditions: application to linear-scaling density functional theory.

We present a comparison of methods for treating the electrostatic interactions of finite, isolated systems within periodic boundary conditions (PBCs), within density functional theory (DFT), with particular emphasis on linear-scaling (LS) DFT. Often, PBCs are not physically realistic but are an unavoidable consequence of the choice of basis set and the efficacy of using Fourier transforms to compute the Hartree potential. In such cases the effects of PBCs on the calculations need to be avoided, so that the results obtained represent the open rather than the periodic boundary. The very large systems encountered in LS-DFT make the demands of the supercell approximation for isolated systems more difficult to manage, and we show cases where the open boundary (infinite cell) result cannot be obtained from extrapolation of calculations from periodic cells of increasing size. We discuss, implement, and test three very different approaches for overcoming or circumventing the effects of PBCs: truncation of the Coulomb interaction combined with padding of the simulation cell, approaches based on the minimum image convention, and the explicit use of open boundary conditions (OBCs). We have implemented these approaches in the ONETEP LS-DFT program and applied them to a range of systems, including a polar nanorod and a protein. We compare their accuracy, complexity, and rate of convergence with simulation cell size. We demonstrate that corrective approaches within PBCs can achieve the OBC result more efficiently and accurately than pure OBC approaches.